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Abstract 

We  study  the  dynamics  and  explore  the  controllability  of  a  family  of 
sphere-plate  mechanical  systems.  These  are  nonholonomic  systems  with  a 
five-dimensional  configuration  space  and  three  independent  velocities.  They 
consist  of  a  sphere  rolling  in  contact  with  two  horizontal  plates.  Kinematic 
models  of  sphere-plate  systems  have  played  an  important  role  in  the  control 
systems  literature  addressing  the  kinematics  of  rolling  bodies,  as  well  as 
in  discussions  of  nonholonomic  systems.  However,  kinematic  analysis  falls 
short  of  allowing  one  to  understand  the  dynamic  behavior  of  such  systems. 

In  this  work  we  formulate  and  study  a  dynamic  model  for  a  class  of  sphere- 
plate  systems  in  order  to  answer  the  question:  “is  it  possible  to  impart  a 
net  angular  momentum  to  a  sphere  which  rolls  without  slipping  between 
two  plates,  given  that  the  position  of  the  top  plate  is  subject  to  exogenous 
forces?” 

*This  work  was  funded  by  ARO  grant  DAAG  55  97  1  0114. 

Uhls  work  was  done  while  the  author  was  with  the  Division  of  Engineering  and  Applied  Sciences  at 
Harvard  University,  Cambridge,  MA. 
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Figure  1:  A  Sphere-Plate  System. 


1  Introduction 

In  his  influential  book  on  analytical  dynamics,  E.T.  Whittaker  [1]  brought  together  many  of 
the  significant  results  in  classical  mechanics,  up  to  that  time.  He  included  a  rich  collection  of 
problems  in  rigid  body  dynamics,  some  involving  nonholonomic  systems  and  rolling  motion. 
More  recently,  with  applications  of  mathematics  branching  out  to  the  areas  of  robotics  and 
object  manipulation,  such  systems  have  come  to  the  fore  yet  again,  this  time  with  emphasis 
on  control  theory.  In  this  work  we  make  a  detailed  study  of  what  is  perhaps  one  of  the 
simplest  nonholonomic  mechanical  systems,  called  the  sphere-plate  system.  The  system 
under  consideration  consists  of  a  sphere  rolling  between  two  horizontal  plates.  The  bottom 
plate  is  regarded  as  being  fixed,  while  the  top  plate  is  movable  in  the  horizontal  plane  (see 
Figure  1).  Our  purpose  is  twofold.  We  want  to  first  formulate  a  dynamic  model  for  a  family 
of  sphere-plate  systems  and  then  to  use  that  model  in  the  context  of  control  theory,  in  order 
to  explore  the  controllability  properties  of  these  systems. 

It  is  intuitively  known  —  and  can  be  proven  mathematically  —  that  a  sphere  can  be 
arbitrarily  re-positioned  and  re-oriented  on  a  plane  by  rolling  (see  for  example  [2],  [3],  [4]). 
In  this  paper  we  investigate  the  problem  of  spinning-up  the  sphere  by  moving  the  top  plate. 
By  the  term  “spin”  we  understand  the  rotational  velocity  of  the  sphere  about  the  vertical 
axis  through  its  center.  The  dynamic  model  will  allow  us  to  define  a  “zero-translational- 
velocity  submanifold”  as  the  set  of  system  states  that  correspond  to  a  spinning  sphere  with 
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the  top  plate  being  stationary.  This  idea  will  be  formalized  in  Sec.  2.3.  What  we  find  is  that  a 
fundamental  relationship  exists  between  the  inertial  symmetry  of  the  sphere  and  the  problem 
of  achieving  a  desired  spin.  In  that  respect,  the  class  of  sphere-plate  systems  considered  in 
this  work  exhibit  the  same  critical  reliance  on  inertial  asymmetry  and  nonholonomic  effects 
as  the  recently  popular  “rattleback  top”  (see  [5],  [6],  [7]).  Specifically,  we  show  that  the  only 
spheres  to  which  one  can  impart  a  steady  spin  are  those  whose  geometric  centers  and  centers 
of  mass  do  not  coincide.  Although  we  focus  on  a  specific  nonholonomic  system,  the  study  of 
such  questions  can  contribute  towards  understanding  more  general  spin- up  problems  [8]. 

The  kinematics  of  sphere-plate  systems  are  a  prototype  for  more  general  nonholonomic 
systems  and  as  such  they  have  been  used  to  demonstrate  key  ideas  in  nonlinear  control. 
In  [9],  [10],  optimal  control  problems  were  formulated  for  a  simplified  kinematic  model, 
sometimes  referred  to  as  the  “nonholonomic  integrator”  (see  Appendix).  More  recently,  [11] 
used  the  same  model  to  illustrate  new  ideas  in  pattern  generation  and  approximate  inversion. 
Kinematic  models  for  sphere-plate  systems  have  also  received  attention  in  previous  works  on 
the  control  problem  of  re-positioning  and  re-orienting  rigid  bodies  under  rolling  constraint 
([2],  [12],  [13]  and  others).  The  ideas  of  reachability  and  Lie  algebras  play  a  natural  role  in 
that  setting.  In  [14],  [2]  and  others,  algorithms  were  presented  for  deciding  the  existence  of 
admissible  paths  between  contact  configurations  of  two  rolling  bodies  and  for  finding  such 
paths.  Those  algorithms  were  in  turn  applied  in  the  areas  of  robotics  and  multi- fingered 
manipulation  ([15],  [16],  [17],  [18]  and  others). 

To  advance  our  understanding  of  sphere-plate  systems  beyond  what  is  afforded  by  kine¬ 
matics,  we  will  draw  on  [1]  which  discusses  a  general  approach  to  the  equations  of  motion  for 
nonholonomic  systems  (originally  due  to  Hamel),  as  well  as  a  number  of  interesting  examples 
relevant  to  this  work.  Our  approach  is  consistent  with  nonholonomic  mechanics  where  the 
equations  of  motion  are  extrema  of  some  energy  functional,  allowing  arbitrary  variations  on 
the  coordinates  and  using  Lagrange  multipliers  to  effect  the  constraints.  This  is  in  contrast 
to  the  vakonomic  approach  [19]  where  the  energy  function  is  first  restricted  to  a  submanifold 
defined  by  the  constraints  and  variations  are  allowed  only  on  that  submanifold. 
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In  engineering,  one  often  wants  to  consider  the  mechanical  system  not  just  from  the 
point  of  view  of  analytical  mechanics  but  also  from  the  point  of  view  of  control  systems.  In 
analytical  mechanics  the  equations  of  motion  for  a  dynamical  system  are  obtained  with  em¬ 
phasis  on  describing  the  behavior  of  the  mechanical  system.  In  the  control  systems  approach 
one  seeks  instead  to  identify  exogenous  input  signals,  obtain  input/output  formulations  and 
prescribe  desired  system  responses  [20].  Both  of  these  points  of  view  are  relevant  to  the  goals 
of  this  work  and  will  be  used  as  appropriate  in  our  discussion. 

In  Sec.  2  we  derive  the  equations  of  motion  for  a  class  of  sphere-plate  systems.  We 
show  that  asymmetry  plays  a  critical  role  with  respect  to  the  controllability  and  existence 
of  integral  invariants  for  the  sphere-plate  system,  viewed  as  a  control  system.  By  the  term 
“integral  invariant”  we  understand  a  function  that  depends  only  on  the  state  of  the  control 
system  and  has  constant  value  along  trajectories  regardless  of  the  choice  of  control  inputs.  In 
Sec.  3  we  present  simulation  results  that  demonstrate  some  control  strategies  for  achieving 
spin.  The  choice  of  control  inputs  that  produced  spin  was  guided  by  physical  intuition  about 
the  mechanical  system.  It  would  be  desirable  to  compute  inputs  that  achieved  a  desired  spin 
and  also  satisfied  some  type  of  optimality  condition,  such  as  minimum  energy  or  minimum 
time,  to  achieve  that  spin.  That  problem  remains  open. 


2  A  Class  of  Asymmetric  Sphere-Plate  Systems 

Consider  a  sphere  of  unit  radius  (Figure  2)  that  rolls  without  slipping  between  two  plates, 
both  plates  being  horizontal  with  respect  to  some  spacehxed  inertial  frame.  It  is  assumed 
that  both  sphere-plate  contacts  (top  and  bottom)  are  maintained  at  all  times.  The  bottom 
plate  is  fixed  while  the  top  plate  is  allowed  to  move  horizontally,  acted  on  by  external  forces 
or  by  a  potential  held.  We  will  ignore  gravity. 

The  five-dimensional  configuration  space  for  the  sphere-plate  system  is  C  =  M2  x  SO( 3). 
The  phase  space  is  three-dimensional  due  to  the  rank-2  rolling  constraint  that  is  imposed. 
We  choose  coordinates  on  the  space  C  as  follows:  Let  T  be  a  bodyhxed  inertial  frame  whose 
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Figure  2:  A  Class  of  Sphere-Plate  Systems. 

origin  is  fixed  at  the  center  of  the  sphere.  The  matrix  0  G  SO  (3)  will  describe  the  orientation 
of  the  sphere.  The  columns  of  0  are  the  spacehxed  coordinates  of  the  unit  vectors  of  T .  The 
vector  x  =  [x\,  ]T  G  M2  specifies  the  horizontal  position  of  the  center  of  the  sphere.  Unless 
otherwise  noted,  quantities  will  be  expressed  in  a  spacehxed  coordinate  frame  whose  z-axis 
is  normal  to  the  plates.  A  potential  held  U(x)  acts  on  the  top  plate,  causing  it  to  move 
horizontally.  We  will  take  U  to  be  differentiable.  The  sphere-plate  system  can  be  considered 
an  autonomous  dynamical  system  with  the  top  plate  being  influenced  by  the  potential  U(x). 
From  a  control  system  point  of  view,  we  will  take  a  pair  of  external  forces  U\,U2  acting 
horizontally  on  the  top  plate  to  be  the  exogenous  inputs,  with  the  potential  held  U  being 
identically  zero. 

In  the  bodyhxed  frame  T  the  sphere’s  center  of  mass  has  coordinates  d  G  M3,  with 
\d\  =  r  <  1.  The  mass  of  the  sphere  is  ms  >  0.  The  rotational  inertia  measured  about  any 
axis  through  the  center  of  mass  is  ams  with  a  >  0.  The  mass  of  the  top  plate  is  mp  >  0. 
Under  these  assumptions,  the  rotational  inertia  of  the  sphere,  expressed  in  any  frame  located 
at  its  center  of  mass,  is  given  by  the  3x3  matrix  amsI.  We  take  e«  to  be  the  ith  standard 
basis  vector  in  Mn  and  dehne  the  quantities: 

dt  =  Od-e3 

a  (!) 

dfr  —  Qd  T  e3 

which  represent  the  (spacehxed)  vectors  from  the  top  and  bottom  contacts  respectively  to 
the  sphere’s  center  of  mass.  In  the  following,  numerical  subscripts  will  be  used  to  indicate 
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elements  of  vectors,  unless  otherwise  noted.  For  example,  x2  =  {x,e2)- 

Armed  with  the  above  definitions,  we  will  first  obtain  the  equations  of  motion  for  the 
sphere-plate  system  assuming  the  mass  of  the  top  plate  is  mp  =  0.  We  will  subsequently 
augment  our  model  to  include  the  inertial  effects  of  the  top  plate. 


2.1  Case  I:  Top  Plate  with  zero  Mass  (rnp  =  0) 

Let  uo  G  M3  be  the  angular  velocity  of  the  sphere,  expressed  in  spacehxed  coordinates.  The 
position  and  orientation  of  the  sphere  evolve  according  to 

=  S(uj)e3  (2) 

0 

0  =  S(u)Q  (3) 


where  S  :  M3  — >  so(3)  is  a  transformation  that  takes  vectors  in  M3  into  3x3  skew-symmetric 
matrices  1 .  Note  that  Eq.  2  expresses  the  rolling  constraint  which  dictates  that  the  instan¬ 
taneous  velocity  of  the  bottom  contact  should  be  zero. 

Let  p  G  M3  be  the  spacehxed  position  of  the  center  of  mass  of  the  sphere.  If  KE  is  the 
sphere’s  kinetic  energy,  then  the  equations  of  motion  for  the  sphere  are  [1]: 


d  dKE 
dt  dp 
d  dKE 
dt  duj 


+  7  x 


dKE 

dp 


+  7  x 


dKE 

du> 


dU(x) 

dp 


T-S 


(  dU(x )  \ 
\  dp  ) 


dt 


(4) 

(5) 


where  7  G  l3  is  the  vector  of  angular  velocities  of  the  inertial  frame  in  which  we  chose 
to  express  matters.  In  our  case,  7  =  0  because  the  spacehxed  frame  does  not  rotate.  The 
vectors  F,T  €  M3  are  the  external  forces  and  torques  (including  those  forces  necessary  to 
1For  x,  y  €  M3,  x  x  y  =  S(x)  ■  y  with 

0  —£3  x2 

s(x)  =  £3  0  -£1 

—£2  x\  0 
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enforce  the  nonholonomic  constraints)  acting  on  the  center  of  mass  of  the  sphere.  Note  that 
we  have  chosen  to  superimpose  the  effects  of  the  potential  U  together  with  those  of  the 
exogenous  forces  u  (implicitly  included  in  F  and  T)  that  are  acting  on  the  top  plate.  This 
is  done  for  convenience  since  we  can  subsequently  set  either  U  or  «  to  zero,  thus  obtaining 
the  equations  of  motion  for  the  control  system  or  the  autonomous  system  respectively. 

Assuming  the  mass  of  the  top  plate  is  zero,  the  kinetic  energy  for  the  system  is 


KE  =  ~^(pTp  +  awTw) 

Zi 


(6) 


with  the  linear  velocity  of  the  sphere’s  center  of  mass  given  by 


P  = 


x 

0 


+  S(u)Qd  =  S(u)db 


(7) 


Substituting  for  the  kinetic  energy  and  linear  velocity  (Eq.  6,  7)  into  the  equations  of 
motion  (Eq.  4,5)  we  obtain  a  set  of  equations  that  involve  the  evolution  of  the  angular 
velocities: 


ms  {S(u)db  +  S(u)2Od)  =  F-^ 


(8) 


amsuj  =  T  -  S  dt 

The  external  forces  acting  on  the  sphere  do  so  only  through  the  contacts  with  the  plates. 
Let  ttjCGl3  be  the  external  forces  acting  at  the  top  and  bottom  contacts  respectively.  The 
vector  c  corresponds  to  the  forces  that  are  necessary  to  enforce  the  rolling  constraint.  If  we 
express  u  and  c  in  spacehxed  coordinates,  we  can  write: 


F  =  c  +  u,  T  =  S(u)dt  +  S(c)db 


(9) 


Because  u3  and  c3  are  co-linear  and  the  vectors  dtl  dh,  have  the  same  projection  onto  the 
horizontal  plane,  we  can  combine  the  effects  of  U3  with  c3  and  regard 

u  =  [ui  U2  0]T  (10) 

to  be  the  external  force  applied  to  the  top  plate. 
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We  combine  Eq.  8,  9,  into  a  system  of  six  equations 


(11) 


ms  (S(Co)db  +  S(uo)20d)  =  u  +  c  - 

amsCo  =  s(u-^jjfydt4rS(c)db 
that  must  be  solved  for  the  evolution  of  to.  Eliminating  the  constraint  forces  c,  we  obtain: 


(al  -  S(db)2)  Co  =  ^-S(e3)  hi  “  j  “  S(db)S(oo)2<dd  (12) 

If  we  choose  to  express  the  sphere’s  orientation  using  roll-pitch- yaw  angles  then  for  6  G  K3 
we  can  write: 

0  =  Rotz{93)  Roty(02)  Rotx(Qi)  (13) 

Without  loss  of  generality,  we  can  take  the  center  of  mass  to  be  located  along  the  cc-axis  of 
the  bodyhxed  frame  and  write  d  =  re i,  where  r  =  \d\.  In  that  case,  we  observe  that  for 
a  >  0, 

det(al  —  S(db )2)  =  a(l  +  a  +  r2  —  2r  sin  d2)2  >  0  (14) 

This  determinant  is  positive  as  long  as  r  <  1,  so  that  the  center  of  mass  is  located  inside  the 
sphere.  In  particular,  the  matrix  ( al  —  S(db)2)  is  symmetric,  positive  definite  for  a  >  0  and 
r<l,  with  eigenvalues 

A1;  A2  =  1  +  a  +  r2  —  2r  sin  92,  A3  =  a.  (15) 


We  summarize  the  equations  of  motion  for  the  rolling  sphere: 


x 

0 

Co 


U02 

-Ml 

S(uo)0 

(al  -  Sfd,)2)-1  (+S(ea)  (u  -  f )  -  S(db)S(^Bd  ) 


(16) 


Alternatively,  using  the  orientation  angles  0,  the  equation  for  the  evolution  of  0  (second  of 
Eq.  16)  can  be  replaced  by: 


6  =  K(0)u 


(17) 


with 


sec  02  cos  #3  sec  62  sin  0:i  0 

m  =  —  sin  6*3  cos  93  0  (18) 

cos  93  tan  ()•>  sin  0:>  tan  d2  1 

2.2  The  Sphere-Plate  Equations 

We  will  now  modify  the  model  presented  in  the  previous  Section,  to  include  the  effects  of  a 
top  plate  with  mass  mp  >  0.  For  this  purpose,  we  consider  the  interconnection  of  the  plate 
and  sphere  systems.  We  note  that  the  top  plate  has  a  translational  velocity  which  is  twice 
that  of  the  center  of  the  sphere.  If  us  G  R2  is  the  horizontal  force  applied  to  the  sphere  by 
the  top  plate  and  up  =  [u\  U2]T  G  R2  is  the  external  force  acting  horizontally  on  top  plate, 
then  the  dynamics  of  the  rolling  sphere  (Eq.  16)  combined  with  the  dynamics  of  the  top 
plate 

2  mpx  =  up  —  us  (19) 

lead  to  the  equations  of  motion  for  the  sphere-plate  system: 

U>2 

X  = 

-LVi 

0  =  S(u)Q 

d,  =  J-'(-S(dt)S(u,fed+  £S(es)  («  -  f )) 

with  u  =  [tii  0]T  and 


We  will  call  the  above  set  of  equations  (Eq.  20,  21)  the  “sphere-plate  equations”.  In  the 
special  case  where  U  =  0  we  will  refer  to  them  as  the  “sphere-plate  control  equations”.  We 
remark  that  J  has  units  of  length  squared.  It  can  be  checked  that  the  matrix  J  is  symmetric, 
positive  definite  for  a  >  0. 
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2.3  The  Control-System  Viewpoint 

We  are  now  ready  to  explore  the  dynamic  model  given  by  the  sphere-plate  control  equations 
(Eq.  20)  and  address  the  question  of  whether  or  not  it  is  possible  to  force  the  sphere  to 
spin  without  translational  motion  by  some  choice  of  control  inputs  U\,U2-  In  the  following 
we  have  in  mind  that  there  are  no  potential  energy  terms  (i.e.  U  =  0)  so  that  the  control 
inputs  are  the  only  external  forces  acting  on  the  top  plate.  We  will  consider  the  special  case 
mp  =  0  for  the  top  plate.  For  a  top  plate  with  mass  rnp  >  0  the  discussion  is  largely  similar 
but  with  more  cumbersome  arithmetic.  That  case  will  be  omitted  here.  However,  the  case 
mp  >  0  will  be  addressed  in  our  simulation  of  the  sphere-plate  system  presented  in  Sec.  3. 

If  we  define  the  system  state  to  be  the  vector 

x 

d 
uo 

with  y  G  R2  x  TSO( 3),  then  we  can  re-write  the  sphere-plate  control  equations  in  the  form 

X  =  fix)  +  Y,9i(x)ui  (22) 

1=1 

where  /  G  R8  is  the  drift  and  gi ,  ,g2  G  R8  the  control  vector  fields. 

We  are  interested  in  knowing  whether  or  not  there  exist  control  inputs  that  can  steer  the 
system  from  any  initial  state  to  any  other  state,  i.e.  connect  a  pair  of  states  by  a  trajec¬ 
tory  in  the  state-space.  For  this  purpose,  we  consider  the  reachable  set  associated  with  the 
sphere-plate  control  system.  For  the  sake  of  completeness,  we  include  some  of  the  relevant 
definitions.  Detailed  discussions  of  the  following  concepts  can  be  found  in  [3],  [4],  [21]  and 
others. 


A 

X  = 


Definition  1  The  reachable  set  associated  with  a  control  system  is  the  set  of  states  that  can 
be  reached  by  the  system  starting  from  an  arbitrary  initial  condition  and  using  appropriate 
control  inputs. 
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Definition  2  A  control  system 


X  =  f(x)  +  J29i(x)ui  (23) 

i=  1 

with  x  in  an  smooth  n-dimensional  manifold  X  and  u  G  Mm  is  controllable  if  for  any  two 
states  XojXi?  there  exist  a  finite  time  T  and  control  inputs  tq(-)  defined  on  [0,T]  so  that 

ft=T  /  m  \ 

Xo  +  jf_o  yf(x)  +  J29i(x)ui)  dt  =  Xi 

Definition  3  The  controllability  Lie  algebra  associated  with  the  control  system  of  Eq.  23  is 
the  Lie  algebra  generated  by  the  drift  and  control  vector  fields:  {f,gi,  ■  ■  ■  ,gm}LA- 

The  sphere-plate  system  evolves  on  a  subset  T>  of  the  ten-dimensional  tangent  bundle 
TC,  C  being  the  five-dimensional  configuration  space.  Because  the  rolling  constraint  (Eq.  2) 
is  of  rank  2  everywhere,  T>  =  M2  x  TSO( 3)  is  an  eight- dimensional  submanifold  of  TC.  In  T> 
we  define  the  following  two  submanifolds: 

Definition  4  The  zero-translational-velocity  submanifold  Z  is  the  manifold 
Z  =  {x  e  V,  (x,  e6)  =  (x,  e7)  =  0}  =  M2  x  SO( 3)  x  S'1 
Definition  5  The  configuration  submanifold  K.  is  the  manifold 

!C  =  {xeZ,(x,es)  =  0}  =  R2xSO(3) 

The  set  Z  contains  state  values  for  which  the  top  plate  is  stationary.  It  is  a  six- 
dimensional  submanifold  of  T>  because  the  constraint  x  =  [0  0]  has  rank  2  everywhere 
on  T>.  Similarly,  /C  is  defined  by  imposing  a  rank-1  constraint  on  Z.  Therefore  K,  is  a  five¬ 
dimensional  submanifold  of  Z  and  is  isomorphic  to  the  configuration  space  C.  We  observe 
that  /C  C  Z  C  T>.  It  is  known  that  1C  is  reachable.  In  other  words,  the  sphere  can  be  arbi¬ 
trarily  re-positioned  and  re-oriented  by  appropriate  top  plate  motions:  this  fact  corresponds 
to  the  kinematic  controllability  of  the  sphere-plate  system.  In  the  following  we  show  that 
the  sphere-plate  system  is  controllable  in  V  if  and  only  if  r  >  0. 
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Because  the  reachable  set  associated  with  a  control  system  is  invariant  under  feedback, 
we  can  apply  a  transformation  to  the  controls,  of  the  type 

u  (->•  ip(x)  +  u  (24) 

ip  :  V  — *  M2  continuous 

in  order  to  cancel  the  drift  terms  in  the  evolution  of  without  altering  the  controllability 

properties  of  the  system.  This  cancellation  of  the  drift  terms  for  ooi,u2,  is  possible  because  if 
r  <  1,  the  3x3  matrix  J  that  enters  the  sphere-plate  dynamics  (last  of  Eq.  20),  is  positive 
definite.  Thus,  its  upper-left  2x2  sub-matrix  is  invertible,  allowing  us  to  solve  for  the  inputs 
that  force 

(hi  =  Ui,  002  =  U2 


Without  loss  of  generality,  we  choose  d  =  re i  for  the  bodyhxed  location  of  the  center  of 
mass  relative  to  the  center  of  the  sphere.  The  resulting  expressions  for  the  drift  and  control 
vector  fields  are: 


U>2 

0 

0 

-uq 

0 

0 

K(6)oo 

,  9i  = 

0 

,  92  = 

0 

0 

1 

0 

0 

0 

1 

"f 

.  u9l  _ 

.  U92  _ 

(25) 


where 


vf 

A 

2  ujtQuj 

r  cos  d2 — - — 

0 

vgi 

A 

2 r  cos  9 2  cos  93  (1  —  r  sin  02) 

"92 

A 

2 r  cos  02  sin  9^(1  —  r  sin  d2) 

8 

A 

2a  +  r2  +  r2  cos(2d2) 

(26) 
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and  Q  is  a  rank-2  3x3  matrix: 

2  sin  03  cos  02  cos  #3  —  cos  02  cos(203)  —  sin  02  sin  03 

Q=  —  cos  02  cos (203 )  — 2  sin  03  cos02  cos  03  sin  02  cos  03  •  (27) 

—  sin  02  sin  03  sin  02  cos  03  0 

We  note  that  the  coordinates  x  and  03  are  ignorable,  i.e.  they  do  not  appear  on  the  right- 
hand  side  of  the  equations  of  motion.  This  is  because  we  chose  d  =  re  1  and  because  we 
arranged  matters  so  that  the  rotational  inertia  of  the  sphere  has  the  same  value  about  any 
axis  through  the  sphere’s  center  of  mass.  As  a  result,  rotation  about  the  axis  connecting  the 
geometric  center  and  the  center  of  mass  —  precisely  what  03  measures  —  leaves  the  dynamics 
invariant.  By  ignoring  x  and  0 1,  one  can  consider  a  reduced,  five-dimensional  version  of  the 
sphere-plate  control  system,  with  state  (02,03,u;).  We  will  show  controllability  for  the  full 
eight-dimensional  system.  We  remark  that  the  evolution  equation  for  the  angular  velocities 
(last  of  Eq.  20)  appears  to  be  more  complicated  than  the  rest  of  the  equations  of  motion. 
This  seems  to  be  a  consequence  of  choosing  to  place  the  sphere’s  center  of  mass  away  from 
its  geometric  center.  The  evolution  equation  for  Co  is  simplified  greatly  if  r  =  0  instead, 
however  in  that  case  one  can  check  that  u;3  =  0  and  the  system  is  uncontrollable. 

The  equations  of  motion  reveal  that  the  two  control  vector  fields  commute  and  thus  they 
can  be  thought  of  as  spanning  the  tangent  space  of  a  two-dimensional  submanifold  in  T>. 

Observation  1  Let  /,  <71,772  G  K8  be  the  drift  and  control  vector  fields  corresponding  to 
the  sphere-plate  control  equations  (Eq.  25).  Then,  [<71,  <72]  =  0  and  Intfigi.g-fi  is  a  two- 
dimensional  submanifold  ofV. 

Proof:  The  drift  and  control  vector  fields  corresponding  to  the  dynamic  model  of  the 
sphere-plate  system  are  given  in  Eq.  25.  The  statement  [<71,(72]  =  0  can  be  checked  directly 
from  the  algebraic  expressions  of  <71,(72-  Alternatively,  notice  that  the  vector  fields  /,  g  1 ,  (.]-> 
are  of  the  form: 


AM 

0 

0 

Mw)  = 

/2(0,(n)  ,  0i  (0)  = 

0  ,  02(0)  = 

0 

Js(0,co)\ 

h(0)e 1 

h(0)e2 
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(28) 


where  /i  6  l2,  /2,/3£  M3  and  6  is  the  vector  of  the  sphere’s  orientation  angles.  The  3x3 
matrix  h(6)  enters  only  in  the  u  equation  but  does  not  depend  on  u.  From  this  fact  and 
from  the  definition  of  the  Lie  bracket,  we  conclude  that  \g\ ,  <72]  =  0,  and  that  all  elements  in 
{f,  91,92} la  involving  [gi,g2]  are  identically  zero. 

The  above  statement  implies  that  at  least  locally,  there  exists  a  coordinate  change  under 
which  both  gi  and  g2  are  constant,  therefore  there  exists  a  two-dimensional  integral  manifold 
whose  tangent  space  is  spanned  by  g\,g2.  We  observe  that  gi,g2  are  independent  and  nonzero. 
Because  SO( 3)  is  parallelizable,  we  have  TSO{ 3)  =  SO{ 3)  x  M3.  Therefore,  Int(gi,g2 )  is  a 
two-dimensional  submanifold  of  SO( 3)  x  M3.  ■ 


Proposition  1  Consider  a  sphere-plate  control  system,  with  a  unit  sphere  of  mass  ms  and 
a  top  plate  of  mass  mp.  Let  ams  be  the  sphere’s  rotational  inertia  about  any  axis  through  its 
center  of  mass.  Furthermore,  let  del3  be  the  location  of  the  center  of  mass  relative  to  the 
center  of  the  sphere,  with  \d\  =  r  <  1. 

For  x  G  R2,  0  G  SO(3),  w£l3  and  u  =  [u\  u2  0]T,  the  set  of  equations 


x 

0 

to 


U)2 

-W 1 


S(u)Q 

J~l  (-S(Qd  +  e3)S(uj)20d  +  —S(e3)u 


describes  the  motion  of  the  sphere-plate  system.  This  two-input  system  is  controllable  in  the 
eight- dimensional  space  V  if  and  only  if  r  >  0,  i.e.  if  and  only  if  the  sphere’s  center  of  mass 
does  not  coincide  with  its  geometric  center. 


Proof:  If  r  =  0  then  the  center  of  mass  is  located  on  the  segment  joining  the  top 

and  bottom  sphere-plate  contacts.  The  torque  applied  to  the  sphere  about  the  spin  axis 
is  identically  zero.  Therefore  the  angular  momentum  of  the  sphere  about  the  spin  axis  is 
conserved  and  is  an  integral  invariant  for  the  system.  We  conclude  that  the  system  cannot 
be  controllable. 
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If  r  >  0,  let  f(6,u)  G  R8  be  the  drift  and  gi(9),g2(9)  G  R8  be  the  control  vector  fields 
associated  with  the  sphere-plate  control  equations,  after  we  have  used  feedback  to  eliminate 
the  drift  terms  in  the  equations  for  chi,  uj2  (see  Eq.25  for  the  case  of  mp  =  0).  If  the  manifold 
V  was  not  reachable,  every  8x8  matrix  P  composed  of  elements  from  {/,  {J\  ■  g2}hA  would 
have  a  determinant  vanishing  on  open  sets.  Because  det(P)  is  an  analytic  function  of  the 
state,  it  cannot  vanish  on  open  sets  without  being  identically  zero.  The  following  set  of 
elements  of  the  taken  from  {f,gi,g2}LA- 
gi->g2i 

[f,9ib  [f,92], 

[f,[f>9i]\,  [/,[/,  02]], 

[/,[/,  [f,9i]]],  [/,[/,[/,  92}}}, 

form  an  8  x  8  matrix  P  whose  determinant  is  non-zero  for  generic  choices  of  system  parameters 
and  state  values.  We  conclude  that  det(P)  can  only  vanish  on  isolated  points,  therefore  the 
above  eight  elements  of  {f,gi,g2}LA  span  almost  everywhere  in  T>.  This  fact  combined 
with  Chow’s  theorem,  tell  us  that  the  Frobenius  manifold  spanned  by  {f ,  gi,  g2} la  is  locally 
diffeomorphic  to  R8  which  precludes  the  existence  of  (non-trivial)  integral  invariants  and 
shows  that  T>  is  reachable.  Reachability  implies  that  there  exists  a  choice  of  control  inputs 
to  steer  the  system  from  any  initial  state  to  any  other  state  in  T>,  thus  the  sphere-plate 
system  is  controllable.  ■ 

The  existence  of  the  integral  invariant  for  r  =  0  implies  that  the  set  of  states  Z  —  /C, 
corresponding  to  non-zero  spin,  is  not  reachable  from  /C.  For  r  >  0,  the  reachable  set  is 
the  entire  manifold  T>.  There  does  not  appear  to  be  a  compact  way  to  characterize  the 
elements  of  {/,  <71,  g2}iA  and  because  their  symbolic  expressions  take  considerable  space  we 
will  not  include  them  here.  Note  that  we  were  able  to  show  controllability  for  the  system 
using  elements  taken  exclusively  from  the  adf- chain  of  {/,  (j\ .  g2} la  (i.e.  elements  of  the  form 
adkgi  where  adfg  =  [/,  g}).  This  brings  up  the  question  of  whether  or  not  the  system  can  be 
linearized,  with  the  linearized  version  being  controllable.  Although  the  distribution  spanned 
by  {.9i;  .92}  is  integrable  as  we  have  shown  (an  immediate  consequence  of  [<71,  <72]  =  0);  it 
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can  be  checked  that  the  distribution  generated  by  5  —  {gi,  g2,  [/,  gf\,  [/,  gfi}  is  not  integrable 
(for  example,  [[/,  <7i],  <72]  ^  span{8}).  Therefore  any  linearization  of  the  sphere-plate  control 
equations  gives  rise  to  an  uncontrollable  system.  Because  the  integrability  test  fails  after 
only  the  first  two  elements  in  the  adj-chain.  namely  g  1 ,  go ,  the  sphere-plate  system  can  be 
thought  of  as  a  second-order  linear  system  coupled  with  a  sixth-order  “genuinely”  nonlinear 
system. 

2.4  Towards  Generality 

Controlling  the  angular  momentum  in  a  sphere-plate  system  can  be  viewed  as  a  special  case 
of  a  more  general  spin- up  problem  [8].  Consider  the  n-dinrensional  control  system: 

m 

X  =  (j)(x)  =  f(x)  +  9i(x)ui  (29) 

1 

with  the  state  evolving  on  a  manifold  X  and  ( x,<f>{x ))  G  TX.  We  will  call  this  the  “base” 
control  system  and  use  it  to  define  a  “derived”  control  system  on  T2X  by 

m 

*  =  /(*)  + £*(*)*  (30) 

1 

with  state  (x,x)  G  TX.  We  could  then  pose  the  question:  “given  that  the  base  system  (Eq. 
29)  is  controllable,  under  what  conditions  is  the  derived  system  (Eq.  30)  also  controllable?” 
Towards  answering  this  question  we  can  offer  the  following  facts: 

Observation  2  The  control  vector  fields  of  the  base  system  (Eq.  29)  defined  on  TX,  com¬ 
mute  with  each  other.  In  addition,  if  the  matrix  [g\, ...,  gm]  has  constant  rank  m  then 
Int(gi, ... ,gm )  is  an  m-dimensional  submanifold  ofTX. 

Proof:  To  show  [g^gfi  =  0  for  all  i,j  G  {1,  use  the  definition  of  the  Lie  bracket 

and  the  fact  that  |&  =  0.  If  the  control  fields  are  all  independent  and  of  constant  rank  then 
Frobenius’  theorem  gives  us  the  existence  of  the  integral  manifold.  ■ 

In  general,  once  we  pass  from  the  base  system  to  the  derived  system  we  can  no  longer 
rely  on  the  control  vector  fields  to  span  additional  directions  by  means  of  Lie  bracketing.  If 
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there  are  m  independent  non-zero  control  fields,  the  remaining  2 n  —  m  directions  must  be 
generated  by  bracketing  with  the  drift  f(x).  In  the  case  of  a  base  system  which  is  linear 
time-invariant  (LTI),  the  situation  is  simple: 


Corollary  1  The  LTI  system  x  =  Ax  +  Bu  x  G  Rn,  u  G  is  controllable  if  and  only  if 
the  system  x  =  Ax  +  Bu  is  controllable. 


Proof:  The  proof  follows  by  the  fact  that  the  system  x  =  Ax  +  Bu  is  controllable  if 
and  only  if  the  matrix 

OB  0  AB  0 

B  0  AB  0  A2B  ... 


is  full  rank. 


3  Simulation  Results 

The  evolution  of  the  sphere-plate  system  was  simulated  for  a  sphere  of  unit  radius  with  mass 
ms  =  2.1  kg  and  a  =  0.5.  The  geometric  center  of  the  sphere  was  initially  at  x  —  [0  0]T. 
The  mass  was  located  at  a  distance  r  =  1  m  from  the  center  of  the  sphere,  along  the  .x-axis 
of  the  bodyfixed  frame.  The  mass  of  the  top  plate  was  mp  =  0.5kg.  Two  simulations  were 
performed,  one  to  spin-up  the  sphere  starting  from  rest  and  the  other  to  increase  the  angular 
momentum  of  an  already  spinning  sphere.  In  both  cases,  the  top  plate  was  required  to  be 
stationary  at  the  beginning  and  end  of  the  simulation.  The  sphere-plate  control  equations 
(Eq.  20)  were  used,  with  the  external  forces  applied  to  the  top  plate  being  the  control  inputs 
and  no  potential  field  present. 

3.1  First  Simulation:  Spin-up  of  a  Stationary  System 

The  system  was  initially  at  rest,  with  the  center  of  mass  located  on  the  sphere’s  equator 
i.e.  0  =  /  or  equivalently,  9  =  0.  A  constant  external  force  u  =  [0  9]TIV  was  applied  to 
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Figure  3:  Spacefixed  Evolution  of  Sphere  -  First  Simulation. 

the  top  plate  for  1.3  seconds  and  then  the  top  plate  was  brought  to  a  stop  by  a  high-gain 
proportional  control  that  servoed  on  the  translational  velocity  of  the  sphere.  Figure  3  shows 
the  spacefixed  evolution  of  the  sphere.  The  center  of  mass  is  shown  as  a  small  dark  ball  on 
the  surface  of  the  sphere.  To  avoid  occluding  the  sphere,  the  top  plate  was  not  drawn.  The 
trajectory  of  the  bottom  contact  is  plotted  on  the  x  —  y  plane.  As  the  sphere  rolls  forward, 
the  center  of  mass  follows  the  trajectory  shown  as  a  cross-hatched  curve.  When  the  top  plate 
comes  to  a  stop  at  the  end  of  the  simulation,  the  center  of  mass  continues  to  move  with  zero 
latitudinal  velocity  and  a  constant  longitudinal  velocity,  corresponding  to  a  steady  spin  for 
the  sphere.  Figure  4  shows  the  trajectory  of  the  center  of  mass  of  the  sphere,  in  a  bodyhxed 
coordinate  frame  whose  origin  is  at  the  center  of  the  sphere  and  whose  axes  remain  parallel 
to  those  of  the  spacefixed  frame.  The  spin  angular  velocity  is  shown  in  Figure  5.  Figure  6 
shows  the  time  history  of  the  kinetic  energy  stored  in  the  sphere.  Once  a  constant  spin  was 
achieved  with  the  top  plate  at  rest,  the  control  inputs  were  such  that  the  top  plate  remained 
stationary.  These  forces  can  be  computed  from  the  equations  of  motion,  as  the  solutions  to 
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Figure  6:  Kinetic  Energy  of  Sphere  -  First  Simulation, 
chi  =  ch2  =  0.  They  are: 


Ml  = 

2  (1  +  sin  62) 

—mscu3 - - - cos  u2  cos  t) 3 

Z 

(31) 

u2  = 

m^(1  +  Sin<ycos0,Smf)3 

Z 

(32) 

3.2  Second  Simulation:  Increasing  the  System’s  Angular  Momen¬ 
tum 

A  second  simulation  was  performed  to  show  an  example  of  increasing  the  angular  momentum 
of  an  already  spinning  sphere.  The  inertia  parameters  were  unchanged  from  the  previous 
simulation.  The  initial  conditions  were  0(0)  =  0  and  u>(0)  =  [0  0  2 ]Trad/s. 

The  strategy  used  to  increase  the  energy  associated  with  the  spinning  sphere  was  moti¬ 
vated  by  physical  intuition:  With  the  top  plate  at  rest,  the  control  inputs  of  Eq.  31,32  will 
keep  the  top  plate  stationary,  (i.e.  chi  =  ch2  =  0).  In  that  case,  the  sphere  will  remain  in  the 
same  horizontal  position  while  spinning  with  a  constant  angular  velocity  u> 3.  While  in  that 
situation,  we  expect  to  be  able  to  “pump”  energy  into  the  system  by  slightly  “leading”  the 
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Figure  7:  Spacefixed  Evolution  of  Sphere  -  Second  Simulation, 
angle  of  rotation  ()■.>;. 


Ui 


u2 


2(1  +  sin.  62 ) 
-mau  3 - - - 

2(1  +  sin  62 ) 
~msu  3 - - - 


COS  02  COs(#3  +  e ) 

cos  02  sin  (03  +  e) 


(33) 

(34) 


In  our  simulation,  we  applied  e  =  0.3 rad  starting  at  t  =  1  sec.  At  t  —  3 sec,  the  top  plate 
was  again  brought  to  rest  using  the  same  high-gain  proportional  control  as  in  the  previous 
simulation.  The  spacefixed  trajectory  of  the  center  of  mass  is  shown  in  Figure  7.  The 
trajectory  of  the  bottom  contact  is  plotted  on  the  x  —  y  plane.  Figure  8,  shows  the  motion 
of  the  center  of  mass  in  the  bodyfixed  coordinate  frame  described  in  Sec.  3.1.  The  center 
of  mass  changes  latitude  as  it  moves,  from  its  initial  location  on  the  equator  to  somewhere 
closer  to  the  pole.  We  would  expect  the  spin  to  increase  as  a  result  of  the  decrease  in  the 
effective  radius  of  rotation.  Aside  from  this  however,  the  increased  spin  was  also  due  to  the 
kinetic  energy  that  was  added  into  the  system  by  the  exogenous  forces.  Figures  9,  10  show 
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X  (m) 

Figure  8:  Bodyfixed  Evolution  of  Center  of  Mass  -  Second  Simulation, 
the  time  histories  of  the  spin  and  the  kinetic  energy  respectively. 

4  Conclusions 

We  have  formulated  a  dynamic  model  for  a  family  of  sphere-plate  systems.  These  are  eight¬ 
dimensional  nonholonomic  systems  that  model  a  sphere  rolling  without  slipping  between 
two  parallel  plates.  Motivated  by  some  recent  questions  on  more  general  spin-up  problems, 
we  have  considered  the  control  problem  of  spinning  the  sphere  between  the  two  plates  by 
means  of  exogenous  forces  applied  to  one  of  the  plates.  The  drift  and  control  vector  fields 
associated  with  the  sphere-plate  control  system  generate  a  Lie  algebra  that  shows  the  system 
to  be  genuinely  nonlinear,  in  the  sense  that  it  cannot  be  approximated  by  a  controllable  linear 
system.  The  control  vector  fields  commute,  as  would  be  the  case  with  any  system  of  the  form 
x  —  f(x)  +  £? gi(x)ui,  with  x  in  a  finite-dimensional  manifold  X.  It  would  be  interesting 
to  explore  more  general  situations  of  this  kind  in  order  to  End  out  under  what  conditions 
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controllability  of  the  system  on  TX  translates  to  controllability  of  the  system  on  T2X ,  as 
suggested  in  [8]. 

The  sphere-plate  system  under  consideration  is  controllable  if  and  only  if  the  sphere’s 
center  of  mass  and  geometric  center  do  not  coincide.  In  that  case,  the  system  can  be 
excited  so  as  to  have  a  spinning  motion  by  an  appropriate  choice  of  inputs.  A  sphere-plate 
system  without  this  inertial  asymmetry  has  an  integral  invariant  associated  with  the  angular 
momentum  of  the  sphere.  We  have  presented  two  control  strategies  for  altering  the  angular 
momentum  of  the  sphere,  starting  from  rest  or  from  a  constant  spin.  Although  we  were  able 
to  find  inputs  that  produced  a  steady  spin  for  the  sphere,  the  problem  of  finding  optimal 
controls  (in  the  sense  of  minimizing  the  integral  of  uTu  or  some  other  cost  functional)  is  still 
open.  The  presence  of  drift  and  the  lack  of  a  compact  characterization  for  the  elements  of 
the  controllability  Lie  algebra  make  this  problem  a  good  candidate  for  applying  some  type 
of  learning  algorithm  in  order  to  obtain  useful  and  efficient  control  inputs. 


Appendix 

Kinematic  models  for  the  sphere-plate  system  are  discussed  in  [2], [10]  and  others.  Briefly,  let 
u,  v  G  R  be  the  velocities  of  the  top  plate,  while  (x,  y )  is  the  position  of  the  bottom  contact 
relative  to  some  coordinate  system  fixed  on  the  bottom  plate.  If  0  G  SO( 3)  describes  the 
orientation  of  the  sphere  then  the  kinematics  of  the  sphere-plate  system  are 


x  =  u 
y  =  v 


{ 

\ 

u 

0  =  -s 

V 

V 

0 

/ 

(35) 


e 


where  S  :  M3  — >  so(3)  is  the  transformation  defined  in  §2.1. 

If  z  is  the  angle  of  rotation  of  the  sphere  about  the  axis  through  the  top  and  bottom 
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contacts  then  the  equations 


x  =  u 

y  =  v  (36) 

z  =  xv  —  yu 

represent  a  simplified  kinematic  model  of  the  sphere-plate  system,  sometimes  referred  to  as 
the  “nonholonomic  integrator”  [10].  The  following  result  shows  that  if  we  require  that  the 
top  plate  be  stationary  at  t  =  0  and  t  =  T  >  0  then  it  is  not  possible  to  change  the  sphere’s 
spin,  z. 

Observation  3  Consider  the  kinematic  model  of  the  sphere-plate  system  given  by  Eg.  36. 
If  the  sphere’s  translational  velocity  is  zero  at  t  =  0  and  t  =  T  >  0  then  i(0)  =  z(T). 

Proof:  From  Eq.  36  we  obtain 


z  —  xv  —  yii  —  xy  —  yx.  (37) 

Because  the  sphere’s  translational  velocity  is  zero  at  t  =  0  and  t  =  T,  we  have  u(0)  =  u(T)  = 
0  and  w(0)  =  v(T)  =  0.  Applying  integration  by  parts  to  Eq.  37  gives: 


(xy  —  yx)dt 


fT  fT 

xy |o  —  xy  dt  —  yx |q  +  /  xy  dt  =  0, 
Jo  Jo 


or  i(0) 


z(T). 
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